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Abstract. A balanced V-shape is a polygonal region in the plane con- 
tained in the union of two crossing equal-width strips. It is delimited by 
two pairs of parallel rays that emanate from two points x, y, are contained 
in the strip boundaries, and are mirror-symmetric with respect to the 
line xy. The width of a balanced V-shape is the width of the strips. 
We first present an 0(n 2 logn) time algorithm to compute, given a set 
of n points P, a minimum-width balanced V-shape covering P. We then 
describe a PTAS for computing a (1 + ^-approximation of this V-shape in 
time 0((n/e) logn + (n/e^ 2 ) log 2 (l/e)). A much simpler constant-factor 
approximation algorithm is also described. 

1 Introduction 

Motivation. The problem we consider in this paper was motivated by the following 
curve reconstruction question: One is given a set of points sampled from a curve 
in the plane. The sample is noisy in the sense that the points lie near the curve, 
but not exactly on it. One would like to reconstruct the original curve from this 
data. Clearly one has to make some assumptions about the point set and the 
curve: If the curve is "too wiggly" or the noise is too large, little can be done. 
One approach is to assume that the curve is smooth and the sample points lie 
not too far from it; see [12,13] and references therein. 1 Roughly speaking, one 
can then approximate a stretch of a curve by an elongated rectangle (or strip) 
whose width is determined both by the curvature of the curve and the amount of 
noise. Refining this approximation allows one to reconstruct the location of the 
curve and its normal vector. 

Complications arise when a curve makes a sharp turn, as it does not have 
a well-defined direction near the point of turn. It has been suggested [12, 18] 

* Work on this paper has been supported by grant No. 2006/194 from the U.S.-Israel 
Binational Science Foundation and by NSF Grant CCF-08-30691. Work by Boris 
Aronov has also been supported by NSA MSP Grant H98230-10-1-0210. An extended 
abstract of this paper appeared in the Proceedings of the Algorithms and Data 
Stuctures Symposium (WADS '11) [7]. 

See [6] for a detailed survey of different notions of measuring similarity between 
geometric objects; is there a sensible (and relevant for our purposes) notion of 
closeness between a discrete unordered point set and a curve? 
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that one approach to handle this situation is to replace fitting the set of points 
corresponding to a smooth arc of a curve with a strip by fitting with a wedge-like 
shape that we call a "balanced V-shape;" perhaps one might incorporate it in 
an algorithm such as that of [16]. It is meant to model one thickened turn in a 
piecewise-linear curve; refer to the figure and precise definitions below. 

In this paper, we construct a slower exact algorithm for identifying a V-shape 
that best fits a given set of points in the plane, then a faster constant-factor 
approximation algorithm, and finally a considerably more involved algorithm 
that produces a (1 + ^-approximation, for any positive e. 

The problem we solve is a new representative of a widely studied class of 
problems, namely geometric optimization or fitting questions; see [1,3-5] and 
references therein. Generally, the problem is to find a shape from a given class 
that best fits a given set of points. Classical examples of such problems are 
linear regression in statistics, the computation of the width of a point set (which 
constructs a minimum- width strip covering the set), computing a minimum 
enclosing ball, cylinder, or ellipsoid, a minimum-width spherical or cylindrical 
shell, or a small number of strips of minimum width, covering the point set; 
see [1,10]. 

Previous work most closely related to our problem is that of Glozman, Ke- 
dem, and Shpitalnik [17]. They compute a double-ray center for a planar point 
set S. A double-ray center is a pair of rays emanating from a common apex, 
minimizing the Hausdorff distance between S and the double ray. While the 
shape they consider is not exactly a V-shape, it is similar enough to be used for 
the same purpose. The exact algorithm they present runs in 0(n 3 a(n) log 2 n) 
time, however, in contrast to our near-quadratic-time algorithm; here a(n) is the 
inverse Ackermann function. 

Another paper closely related to our problem is that of Agarwal, Procopiuc, 
and Varadarajan [2]. It concerns the 2-line-center problem studied extensively 
in the past; see the references in [2]. The goal is to cover a given set of points 
by two strips of minimum possible width. One application is fitting two lines 
to a point set. There had been several previously known near-quadratic-time 
exact algorithms for the problem. An 0(nlogn)-time 6-approximation algorithm, 
and an 0(nlogn + ne~ 2 log(l/e) + e -7 / 2 log(l/e))-time (1 + e)-approximation 
algorithm were presented in [2]. A V-shape covering a point set is a special case 
of covering a point set by two strips, so some of the tools from [2] apply to our 
problem as well. 



Problem statement and results. In this 
paper, we focus on the class of polygonal 
regions in the plane that we call balanced 
V-shapes. A balanced V-shape has two 
vertices x and y and is delimited by two 
pairs of parallel rays. One pair of parallel 
rays emanate from x and y on one side 
of the line xy and the other pair of rays 
emanate from x and y on the other side 
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of xy, symmetrically with respect to xy (see the above figure). In particular, a 
balanced V-shape is completely contained in the union of two crossing strips of 
equal width. Its width is the width of the strips. 

Consider a point set P of n points in the plane. We describe, in Section 3, an 
0(n 2 log n) time algorithm that computes a balanced V-shape with minimum 
width covering P. 

Our algorithm actually identifies a particular type of V-shapes that wc 
call "canonical" (see below for definitions) and enumerates all minimum-width 
canonical V-shapes covering P; as some degenerate n-point sets have 0(n 2 ) such 
V-shapes (see Section 4), this approach will probably not yield a subquadratic 
algorithm. This leaves open the problem of how quickly one can identify just one 
minimum-width V-shape covering P. 

In Section 5, we present an 0(n log n) algorithm that constructs a V-shape 
covering P with width at most 13 times the minimum possible width. In Section 6, 
we show how to construct a (1 + e)-approximation in time 0((n/e) log n + 
(n/e 3 / 2 ) log 2 (l/£)), starting with the 13-approximation obtained earlier. 

2 Reduction to canonical V-shapes 

In the remainder of this paper, for simplicity of presentation and unless noted 
otherwise, we assume that the points of P sue in general position: no three points 
are collinear and no two pairs of points define parallel lines. All algorithms can 
be adapted to degenerate inputs without asymptotic slowdown. 

We will find it convenient to consider a larger class of objects, namely V-shapcs. 
A (not necessarily balanced) V-shape (refer to the figure below) is a polygonal 
region similar to a balanced V-shape ex- 
cept that the widths of its two arms need 
not be the same. More formally, a V- 
shape V is a polygonal region bounded 
by two pairs of parallel rays emanating 
from its two vertices x and y. One pair of 
parallel rays (left rays X\ and Yi) lies on 
the left side of the directed line xy, while 
the other pair (right rays X 2 and Y 2 ) lies 
on its right side. The inner rays Xi emanate from x, while outer rays Yi emanate 
from y. X x U X 2 is the inner boundary of V, while Yi U Y 2 is its outer boundary. 
The left arm of V, L = L(V), is its portion on the left of xy; i.e., it is the 
region bounded by rays X\ and Y\ and segment xy. The width of the left arm, 
width(L(V)), is the distance between Xi and Y\. The right arm and its width 
are defined analogously. The width of V, width(V), is the larger of the widths of 
its two arms. V is contained in the union of two strips Si and S 2 : Si is delimited 
by the lines containing X, and Yi, respectively; we refer to Si and 52 as the left 
and right strip of V, respectively. 

A minimum-width balanced V-shape can be obtained from a minimum-width 
V-shape by widening the narrower arm until the widths of the arms are equal. 
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In the remainder of the paper, the n-point set P is fixed. To avoid trivial 
cases, we assume that n > 5. By the general position assumption, all points of P 
cannot be collinear, nor can P be covered by a V-shape of zero width. We need 
not consider V-shapes with all points in one strip as according to Lemma 1, such 
a V-shape does not have minimum width. 

Lemma 1. In a positive-width minimum-width V-shape V covering P, it is not 
possible that one of the strips already contains P in its entirety. 

Proof. Suppose Si :— Si(V) covers P. Let w > be its width. We argue that 
there is a V-shape covering P of width w/2 + e, for any positive e, so V does not 
have minimum width. Indeed, let £ be the median line of S\. Cut Si by £ into 
two parallel strips of width w/2. They cover P. They do not form a V-shape, but 
they can be approximated arbitrarily closely by a V-shape near P, by placing 
its vertex y sufficiently far to the left of P along £, x at the rightmost point of 
CH(P) n £, and the boundary rays near-parallel to £. □ 

Unless otherwise stated, the only particular V-shapes we will be interested in 
are the ones we call canonical. A V-shape is canonical, if the bounding rays of each 
arm pass through exactly three points of P; more precisely if |XinP| + |YinP| > 3, 
for i = 1,2 (recall that, by our general position assumption, |JQnP|, |5^nP| < 2); 
in addition, we require that each arm of a canonical V-shape covering P is locally 
of minimum width, i.e., neither arm can be narrowed by an infinitesimal motion. 

The reason why we consider only canonical V-shapes is that at least one 
minimum- width V-shape covering P is canonical (see Lemma 2 below), so we can 
confine the search to canonical V-shapes and discard any non-canonical V-shapes 
considered by our algorithm. 

Lemma 2. At least one minimum-width V-shape covering P is canonical. 

Proof. In order to prove this lemma, we first explain why we can assume that 
\Xi n P| + \Yi n P| > 3. We then discuss how the boundary points may be 
positioned on the arms. 

By Lemma 1, in no minimum-width covering V-shape one strip covers P 
completely Hence in the following we disregard this possibility. 

We present a sequence of transformations, starting with a minimum-width 
V-shape covering P, which do not increase its width and end in a canonical 
V-shape. We begin by translating its outer boundary in the direction of first Y\ 
and then Yi to ensure that each of the outer rays contains a point of P; this point 
might be y. Now translating the inner boundary, first in the direction opposite 
to that of X\ , and then that of X 2 , we guarantee that each of X\ , X 2 contains 
a point of P. (By Lemma 1, an outer ray Yi cannot meet its corresponding 
inner ray X^ without meeting a point of P, as we started with a minimum-width 
V-shape.) Now consider an arm (say, left) of the resulting V-shape. We will 
further transform it so that \X\ (~l P| + |Yi fl P| >2. We have so far ensured that 
each of X\ and Y\ contains at least one point. If exactly one point is present on 
each left ray, X\ and Y\ can be rotated so that the width of Si shrinks. This 
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process stops either when Si collapses to a line (in which case it is easy to check 
that Yi contains two points of P and X\ contains at least one) or when three 
points S\, s 2 , s 3 lie in Xi U Yi, say Si and s 2 on one ray and s 3 on the other. In 
the latter case, unless the angles ZssS\S 2 and Zsis 2 s 3 are acute, the rotation 
can proceed, further narrowing S\. Corollary 1 follows from this last condition. 

Repeating the process with the right arm, we arrive at a covering V-shape 
whose width is no larger than that of the original V-shape, with the property that 
(a) it satisfies Corollary 1 if there is no zero-width arm, (b) each bounding ray 
contains a point of P, and (c) |X X n P| + 1 Yi n P| > 3 and \X 2 n P\ + | Y 2 n P\ > 3; 
i.e., the resulting covering V-shape is canonical and as good or better in terms of 
width. Hence, indeed, it is sufficient to examine only canonical V-shapes. □ 

Corollary 1. Let P be a point set and V be a minimum-width canonical V-shape 
covering it. Let si, s 2 , s 3 € X be the points on the boundary of a non- zero-width 
arm of V , with s\ and s 2 on one ray and s 3 on the other. Then the angles 
/s 3 S\S 2 and Zs\s 2 s 3 are acute. 

All canonical minimum- width V-shapes fall into the following three categories: 

both-outer Each outer ray contains exactly two points of P, and each inner 

ray contains at least one. 
inner-outer On one arm of the V-shape, the inner ray contains exactly two 

points of P; on the other arm, the outer ray contains exactly two points. The 

other rays contain at least one point of P. 
both-inner Each inner ray contains exactly two points of P and each outer ray 

contains at least one. 

3 Computing a canonical minimum-width V-shape 

To find a canonical minimum-width V-shape covering P, we will search inde- 
pendently for the best solution for each of the three types identified above and 
output the V-shape that minimizes the width. Let H be the convex hull of P. 

V-shapes of both-outer type. Consider a covering V-shape V with outer rays Y\ , Y 2 
containing edges ei,e 2 of H, respec- 
tively; refer to the figure on the right. 
Let I be the bisector of the angle 
ZY\yY 2 . Notice that V is not mini- 
mal unless its width is given by the 
largest distance from a point in P to 
its closest outer ray. Therefore, we can 
assume that points of P left of I belong 
to the left arm of V and points right of i — to its right arm. 

Thus, given Y\, Y 2 , and i it is sufficient to determine the furthest point 
from Y\ to the left of I and the furthest point from Y 2 to the right of £. The larger 
distance determines the width of V. This can be accomplished by building a data 
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structure D{P) on P that supports the following queries: Given a halfplane h and 
a direction d, return an extreme point of P fl h in direction d. 0(n 2 ) queries are 
sufficient to enumerate all choices of e\ , e 2 and identify the best both-outer-type 
V-shape. D(P) can be constructed in 0(n 2 logn) time and supports logarithmic- 
time queries, resulting in total running time of 0(n 2 logn). 

D(P) is constructed as follows: We build the arrangement A = A{P*) of lines 
dual to points of P. Cells of A correspond to different ways to partition P by 
a line. We construct a directed spanning tree T of the cells of A, starting with 
the bottommost cell and allowing only arcs from a cell / to a cell immediately 
above / and sharing an edge with it; we use Pf C P to denote the convex hull of 
the set of points whose dual lines lie below /. Using T as the history tree, we 
store the convex hull Pf for every face / £ A, in a fully persistent version [14] of 
the semi-dynamic convex hull data structure of [21]. We also preprocess A for 
point location. Given a query (say, upper) half-plane h and direction d, we locate 
the face / of A containing the point dual to the bounding line of h and consult 
the data structure associated with / and storing Pf = P n h to find the extreme 
point of Pf in direction d, all in logarithmic time. 

V-shapes of inner-outer type. In this section, we describe how to find a minimum- 
width canonical V-shape covering P and having exactly one edge of CH(P), say e, 
on its outer boundary; it contains two points of P on the inner bounding ray 
of its other arm. We handle each choice of e independently, in O(nlogn) time, 
yielding overall 0(n 2 logn) running time. 

Having fixed an edge e of CH(P), con- 
sider a (minimum-width canonical) V- £' P1P2 

shape V covering P that has e on its 
boundary. For ease of description, suppose 
Y\ D e, X 2 contains two points Pi,P2 G P, 
while both Y2 and X\ contain at least one 
point of P each, denoted p and q, respec- 
tively; see the figure on the right. 

Let £ be the line containing e, and £' 
be the line containing e' : — P\P2- Set Q := Xi 
Si n P and Q' :=P\ Q. We observe that 

a) Q is the set of points of P at distance y l 
at most dist(q,£) from £; 

b) P1P2 is an edge of CH(Q'); 

c) y 2 is contained in a supporting line 
£ PlP2 of CH(Q') (which must also be a supporting line to CH(P) for V to 
cover P) parallel to £'; this line lies on the same side of £' as Q'; 2 and 

d) width(V) = max(width(S*i),width(S* 2 )) = 
max(dist(<?, t), dist(f , £ PlP2 )). 

Our algorithm enumerates all choices for the point q, in order of decreasing 
distance from £. For the current choice of q, it maintains (the boundary of) 

2 If width(i?(l/)) = 0, we have Q' C £' and Y 2 C £'. 
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CH(Q'), say as an AVL tree, and, for each edge e' of CH(Q')) the distance from 
£' to the furthest point of CH(P) to the right of (i.e., "beyond") £'. Edges 
with distances are stored in a min-heap; the minimum such distance gives the 
minimum width for S2 for the current choice of Si ; the larger of the two widths 
determines the width of the current V-shape. We record the best width of any 
V-shape encountered in the process. 

The algorithm is initialized with the set Q' containing the two points of P 
furthest from £ (the case where Q' contains only one point treated by the both- 
outer case as the zero-width strip S2 can be rotated until it contains one edge of 
CH(P)). A generic step of the algorithm involves moving the current point q from 
Q to Q'. We update the convex hull of Q' by computing the supporting tangents 
from q to the old hull, in O(logn) time. For the two new hull edges e\, e 2 , we 
compute the corresponding supporting lines £ ei ,£ e2 of CH(P), using a suitable 
balanced-tree representation of CH(P), also in logarithmic time. We add the new 
edges with the corresponding widths to the min-heap, after removing from it the 
entries of all the eliminated edges of CH(Q'). The root of the min-heap yields 
the best width for S2 for the current partition {Q, Q'}. The algorithm requires 
presorting points by distance from £ and then a linear number of balanced-search- 
tree and heap operations (since the number of edges inserted is less than 2n and 
each cannot be deleted more than once), for a total running time of 0(n log n) 
for a fixed e, as claimed. 

Working through the entire set P (except for the endpoints of e), in order of 
decreasing distance from £, growing Q' and shrinking Q, we obtain a sequence of 
fewer than n V-shapes which include all the canonical minimum-width V-shapes 
covering P with e on its outer boundary and two other points of P lying on 
the opposite arm's inner boundary (it may include non-canonical V-shapes as 
well, but it is not difficult to check that every combination (e, q, e' , £ e i) examined 
by the algorithm yields a valid V-shape covering P, which is sufficient for our 
purposes). 

To summarize, inner-outer type V-shapes can be handled in total time 
(9(n 2 logn). 

V-shapes of both-inner type. Now a covering V-shape V has points a, b of P on 
its inner ray X\ and points c,d on its inner ray X 2 ; refer to Firgurc 1; points 
a, b, c, d are in convex position, in this counterclockwise order. It is known [20] 
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Fig. 1: Empty wedge defined by ab and cd. 
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that there are at most 0(n 2 ) such wedges W — W(a, b, c, d) determined by 
a quadruple of points a,b,c,d £ P and empty of points of P; note that W 
completely determines V, and, given W, one can construct the corresponding V 
in O(logn) time, so it is sufficient to enumerate all empty wedges W. 

For a pair a, b € P, we compute all pairs c, d, so that W(a, b, c, d) is an empty 
wedge. Let Q(a, b) be the set of all points of P lying to the left of the directed 
line ab. 

Observation 1. W(a, b, c, d), in the above notation, is an empty wedge if and 
only if line cd supports CH(Q) and separates segment ab from Q = Q{a 1 b) (and 
a, b, c, d are in this counterclockwise order). 

Now enumerating all k pairs c, d for a fixed choice of a, b can be done in 
time 0((k + 1) logn), as follows. While handling V-shapes of both-outer type we 
constructed a data structure D(P) which, for a given line (here ab), produces a 
balanced search tree storing the convex hull of the points of P lying to one side 
of the line (here Q = Q(a, b)). Using D(P), we find the point z of Q closest to 
the line ab and traverse the boundary of CH(Q) in both directions from z, to list 
all k edges cd of CH(Q) satisfying the conditions of the above observation. Since 
all such edges are consecutive, it is sufficient to examine k + 2 edges of CH(Q). 
Repeating the procedure for all choices of a, b and recalling that the number of 
empty wedges is at most quadratic, we deduce that the enumeration algorithm 
runs in time 0(n 2 logn). 

4 Maximum number of canonical minimum-width 
V-shapes 

How far is our algorithm from optimality? In Figure 2, starting with the vertex set 
of two congruent regular fc-gons, for a suitably large k, we sketch a construction 
of a set of n points with 0(n 2 ) distinct covering minimum-width V-shapes. The 
idea is that a minimum-width covering V-shape would consist essentially of two 
independently chosen minimum-width strips, each covering one fc-gon. The point 
set is highly degenerate. However, perturbing it slightly yields a point set with 
0(n 2 ) canonical V-shapes with width arbitrarily close to minimum possible. This 
is an indication that any algorithm that explicitly enumerates canonical covering 
V-shapes may have to spend fi(n ) time on this input, thus it is unlikely that 
any algorithm taking our approach can run much faster. On the other hand, for 
this specific input one can encode the 0(n 2 ) optimal V-shapes in 0(n) space. 
This leaves open the possibility that identifying a single minimum-width covering 
V-shape may still be possible in subquadratic worst-case time. 

5 A 13-approximation algorithm 

Given a set of points P, let w be the minimum value such that P can be covered 
by a V-shape of width w. We present an algorithm that computes a V-shape 
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Fig. 2: A sketch of a construction of a set with many minimum- width covering 
V-shapes. Points below and above the two regular fc-gons are added to raise the 
minimum width of any covering V-shape to that of the /c-gons. 

covering P of width at most 13u> in time 0{n log n). For this purpose, we use 
the O(nlogn) time 6-approximation algorithm for the 2-line-center problem 
presented by Agarwal, Procopiuc, and Varadarajan [2]. Recall that the 2-line- 
center problem is the following: Given a set P of n points in IR 2 , cover P by 
two congruent strips of minimum width. We start with the following observation 
which follows from the fact that the union of the two strips of any V-shape 
covering P contains P. 

Observation 2. If w' is the width of two congruent strips of minimum width 
covering P, w' < w. 

Our 13-approximation algorithm proceeds as follows. Use the 6-approximation 
algorithm of [2] to compute two congruent strips of width w" that cover P, with 
w' < w" < Qw'. (It is possible that the two strips just computed are such that a 
V-shape defined by them contains P. In this case we return that V-shape. This 
clearly produces a 6-approximation, due to Observation 2. In the remainder of 
this section, we will assume that this is not the case, in other words, one of the 
two strips has points of P on both sides of it.) Find the median lines £\ and £2 
of the strips. For all points in each strip, project them orthogonally onto £\ and 
£2 respectively (the points in the intersection of the strips are duplicated and 
projected onto both £\ and £2)- Let P' be the resulting set of projected points. 
Compute an exact minimum width V-shape V covering P' (see Section 5.1) in 
O(nlogn) time. The desired approximate V-shape V is obtained by widening V 
by w"/2 in all directions. 

Theorem 3. This algorithm computes a 13-approximation of a minimum-width 
V-shape covering P. 

Proof. Let Vbest be a minimum-width covering V-shape of P, V — a minimum- 
width covering V-shape of P', and V apx — the approximate covering V-shape 



computed by the algorithm. As the points of P have been moved by a distance 
of at most w"/2 to form P' , width(V') < width(Vbcst) + w" . Since V^, px is a 
widened version of V , it contains the points of P. Moreover, width(V r apx ) < 
width(V') + w" < width(Vbcst) + 2w" < w + 12w' < I3w by Observation 2. □ 

Remark. Using the (1 + e)-approximation algorithm of [2] in place of their 
6-approximation algorithm in our procedure, we can attain any approximation 
factor larger than three for the minimum-width V-shape. The running time 
remains O(nlogn), with the constant of proportionality depending on the quality 
of the approximation. We do not discuss this extension further, since we present 
our own (1 + e)-approximation algorithm for the problem in Section 6. 

5.1 Minimum-width V-shape for points on two lines 

We now describe how to compute the minimum-width V-shape V covering a 
given point set P' contained in the union of two lines £\, £2 in the plane. Put 

z :=4n<? 2 . 

Let P[ := P' n £1 and P!, := P'n£ 2 . If P' is covered by a zero- width V-shape, 
which is easy to check, we are done. From now on we assume that this is not the 
case, i.e., that £\, £2 do not already form a V-shape containing P', so £\ separates 
some two points of P'^ and/or £ 2 separates some two points of P[. The convex 
hull CH(P') has three or four vertices. Moreover, by reasoning similar to that 
of Section 3, the outer boundary of V contains two, three, or four vertices of 
CH(P') (in the case where an outer ray is contained in £\ or £ 2 , we consider only 
the extreme points). Before describing how we handle these cases, we need a 
technical lemma. 

Lemma 3. Given a line partitioning P' into P' r , P[ and given their convex hulls 
CH(P^), CR(Pf ), the minimum-width canonical V-shape V of P' containing P' r 
in one strip and P[ in the other can be computed in constant time; some points 
of P' might lie in both strips of V' . 

Proof. The convex hulls CH(P^) and CH(P^) have at most four vertices each. 
It must be the case that the boundary of one arm of V contains an edge e r 
of CH(P^) or an outer common tangent of CH(P^) and CH(P^), and the other 
arm boundary contains an edge ei of CH(P^) or an outer common tangent to 
CH(P^) and CH(P^). There is a constant number of possible pairs of such edges. 
Let S(X, e) be the minimum-width strip covering a set X and parallel to e. 
For each such pair of edges e r , ee, check whether S(P^, e r ) and S{P' t , ei) form a 
V-shape. Return the canonical V-shape of minimum width among all V-shapes 
so generated. □ 

Now we consider the different types of canonical V-shapes covering P' and 
describe how to find a minimum-width V-shape of each type. 

Case 1: An outer bounding ray ofV contains an edge e of CH(P'). Let £ be 
the line containing e. For all points p of P', draw a line £ p through p and parallel 
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to £. Apply Lemma 3 to (the partition induced by) £ p . This can be implemented 
to run in overall time 0(n log n). 

In the remaining cases, each of the outer rays of V contains precisely one 
vertex of CH(P') and each inner ray contains two points of P' . 

Case 2: An inner ray of V lies on l\ or l 2 - Suppose an inner ray of V is 
contained in l\. Draw two lines parallel to l\ and very close to it, one to the left 
of £i, one to the right of l\. Apply Lemma 3 to each of these two lines. 

Case 3: Point z = £\ fl I2 lies between the two arms of V . Draw the two lines 
passing through z and bisecting the angles between £\ and £2- Apply Lemma 3 
to each of these two lines. 

Case 4 : Point z is inside one arm of V . For each pair of consecutive points 
p,q € P' on l\ or on t 2 not separated by z, apply Lemma 3 to the perpendicular 
bisector of the segment pq. 

Now we argue that the last procedure returns the best minimum-width V- 
shapc V' of P' with two points on its outer boundary and z in one of its arms, 
correctly handling case 4 and thereby concluding our description. 

Let x and y be the vertices of V. For ease of presentation, rotate the entire 
picture so that y lies below x; refer to figure 3a. Let S3 be the point of P on Y\, 
Si, S2 be the points on X\, with si closer to x than s 2 . Similarly let S4 be the 
point on Y 2 , s 5 , s 6 be the points on X\, with s 5 closer to x than s 6 . As £\ and 
I2 don't intersect between the two arms of V , si and S5 lie on one line, S2 and 
sq lie on the other line. Let Si and S5 lie on £\, and S2 and sq lie on £2, without 
loss of generality. 

The three points of P' on one arm boundary cannot lie on the same line £j, as 
they form a triangle. Therefore either each line £j contains three boundary points 
belonging to three different boundary rays, or one line contains four boundary 
points from four boundary rays, and the other line contains two boundary points 
from the two inner rays. 

We consider all possible cases: 

a) £1 contains si, S5, and 54, and £2 contains S3, s 2 , and s 6 {£\ contains Si, S3, 
and S5, and £ 2 contains S2, S4, and sq is a symmetric case). 

b) £ 2 contains four boundary points, £\ contains two boundary points. 

c) l 2 contains two boundary points, £\ contains four boundary points. 

In each case, we prove either that the configuration of the points is impossible, 
or that the angles ct\ and ct2 (see figure 3a) between £\ and the inner boundary 
rays are acute. When the angles are acute, the perpendicular bisector of S1S5 
separates the points of P' belonging to each arm of V . Therefore Lemma 3 can 
be applied and our handling of case 4 is justified. 
We consider the cases in turn: 

a) £1 contains Si, s 5 , and S4, and £ 2 contains S3, s 2 , and s 6 (see figure 3a). 
As the angle Zs4s 5 s 6 is acute, its opposite angle ct2 is acute. What is left 
to prove is that ot\ is acute as well. The angle ZS1S2S3 is acute, hence as 
is obtuse, and so is a%. But the angle Zs 5 s 6 S4 is acute, hence £2 does not 
intersect £\ in the right arm of V, so £\ and £2 intersect in the left arm. More 
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(a) Three boundary points on £i,£ 2 - (b) Four boundary points on £ 2 




precisely £\ intersects the segment S3S2. The opposite angle of a.\ is smaller 
than the angle Zs 2 SiS3, which is acute. Therefore so is a.\. 

b) £ 2 contains four boundary points, £\ contains two boundary points. Let £ S2S3 
be the line containing s 2 , S3, £ Se s 4 be the line containing s 6 , 34, and £ xy be the 
line containing x,y (see figure 3b). By Corollary 1, as the angle ZS3S2S1 is 
acute, £ S2 s 3 forms an acute angle a$ with £ xy . Similarly, as the angle Zs^sqs^ 
is acute, £ SeS4 forms an acute angle a± with £ xy . But as £ S3S2 = £ S6S4 , a 3 and 
0:4 are supplementary, a contradiction. 

c) £2 contains two boundary points, £\ contains four boundary points (see 
figure 3c). By Corollary 1, the angles Zs 2 siS3 and /Ls^s^sq are acute, therefore 
their opposite angles a\ and a 2 are acute as well. 

6 A (1 + e) -approximation algorithm 

In this section we describe how to construct, given a point set P and a real 
number e > 0, a V-shape V covering P, with width(V") < (1 + e)u> op t, where 
w opt is the width of a minimum-width V-shape covering P. 

We start by recalling the notion of an anchor pair used in [2]. Given a V- 
shape V covering P, fix one of the strips of V, say Si. We say that a pair of 
points p,q £ P C\ Si is an anchor pair, if dist(p, q) > diam(P n S\)/2. Lemma 3.3 
in [2] describes how to identify at most 11 pairs of points in P, such that, for 
any two-strip cover of P, at least one of the pairs is an anchor pair for one of 
the strips; the algorithm requires 0(n log n) time. As covering by a V-shape is 
a special case of covering by two strips, the definition and the algorithm apply 
here as well. 

We show how to, given a potential anchor pair p, q, construct a (1 + e)- 
approximation of the minimum- width V-shape covering P for which p, q is an 
anchor pair. More precisely, below we prove 

Lemma 4. Given a potential anchor pair p,q e P, we can construct, in time 
0((n/e) logn + (n/e 3 / 2 ) log 2 (l/e)) ; a V-shape covering P, of width at most l + e 
times the minimum width of any V-shape covering P for which p, q is an anchor 
pair. 

Applying this procedure at most 11 times, we obtain our desired approximation 
algorithm: 

Theorem 4. A V-shape covering P and of width at most (1 + e)w pt can be 
constructed in time 0((n/e)logn + (n/s 3 ^ 2 ) log 2 (l/e)). 

We first prove that it is sufficient to consider those V-shapes V with an- 
chor pair p, q, for which the strip containing p, q has one of a small set of 
fixed directions. Setting (3 := sin _1 min{e • width(V)/(6d(p, g)), 1} and 7 := 
^ + sin" 1 ( min{l, width(y)/d(p, g)}), we prove the following 

Lemma 5. Let V-shape V cover P, and let p,q be an anchor pair for Si(V). 
Rotating Si(v) by an angle at most f3 does not increase the width of the V-shape 
by more than a factor of 1 + e/3, and the angle between pq and the direction of 
the rotated strip cannot exceed 7. 
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Proof. Put w := width(y). Let B be the minimum bounding box of PnS\. More 
precisely, it is the shortest rectangle cut out of Si by two lines perpendicular to 
Si and containing PflSi; refer to Figure 4. Let s and t < w be the length (along 
the axis of Si) and width of B, respectively. Let S[ be the minimal parallel strip 




s 



Fig. 4: Rotation by a does not change the width by much. 

containing B n V, whose direction is a < (3 away from that of Si (there are two 
choices for S[, corresponding to rotating clockwise and counterclockwise; only 
one is shown; the argument applies to both cases). Then 

width(S^) < s sin a + t cos a < 2d(p, q) sin a + w 

< w(l + sina) < t(l + e/3), 

w 

since sin a < sin/3 < sw/(6d(p,q)). Now replace Si by S[ to obtain a new V- 
shapc V covering P. Its width is min{width(Si), width(S2)} < (1 + e/3)w, as 
claimed. 

Observe that in the above construction, the angle between pq and the direction 
of S[ cannot exceed 

a + sin _1 (min{l,t/rf(p, q)}) < f3 + sin _1 (min{l, width(V) /d(p, q)}) =7. 

□ 



We conclude that enumerating all V-shapes that contain p, q in their strip Si 
and whose directions are (a) at most 7 away from that of d(p, q) and (b) spaced 
at most (3 apart, would yield a V-shape whose existence is claimed in Lemma 4. 
The number of directions to be tested is at most 0(7/7?) = 0(1 /s). 

Given a candidate anchor pair p, q, the algorithm proceeds by starting with 
the direction pq. Since we need not consider V-shapes whose width is larger 
than the approximate width u> a px computed in Section 5 (this is where the 
13-approximation algorithm is used to bootstrap our (1 + e)-approximation) , we 
replace width(F) by the smaller w ap x/13 in the definition of f3 above and by the 
larger w apx in the definition of 7, thereby erring on the conservative side in each 
case. Having computed (conservative estimates of) /? and 7, we enumerate the 
0(l/e) directions of the form 9i := 9 pq + if}, where 9 pq is the direction of pq and 
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i is an integer ranging from —\"f/(3~\ to [7//?] . It remains to explain how to deal 
with one such direction 9 := 9^. 

Lemma 6. One can compute a canonical V-shape V covering P with one arm 
in given direction 9 and width at most 1 + e/3 times the minimum width of any 
such V-shape, in time 0(nlogn+ (n/e 1 / 2 ) log 2 (l/e)). 

Proof. We use an approach similar to that of the inner-outer case of our exact 
algorithm with a slight twist. 

Let £ be a line in direction 9 supporting CH(P). We again let q be the furthest 
point from I in Q := P n Si and let Q' := P\Q. When q is fixed, the minimum- 
width V-shape is determined by the minimum-width strip S2 covering Q' and not 
"splitting" P, i.e., such that it does not have points of P on both sides of it. It is 
easy to ensure that S2 does not split P by observing that a direction of S2 lying 
between the directions of the common outer tangents to CH(Q) and CH(<3') is 
never useful. Depending on the side where the lines supporting these tangents 
cross, a minimal strip S 2 covering Q' and lying in the range between them either 
crosses Q (and therefore P) or completely covers Q (and therefore P). In the 
former case, Si and S2 do not form a legal V-shape covering P and in the latter 
they form a covering V-shape with one empty strip, which never yields minimum 
width by reasoning as in Lemma 1. 

The width of the resulting V-shape is the maximum of dist (q,£) and (the 
restricted) width(S 2 ). The algorithm proceeds by processing points q in order of 
decreasing distance to £, keeping track of dist(g, I) and a coreset for Q' , which is 
a subset of Q' with the property that its directional width, in every direction, 
is at least 1 — e/3 that of Q' (and, expanding the corresponding minimal strip 
containing the subset by a factor of l + e/3, we get a strip covering Q'). Chan [11], 
in Theorem 3.7 and remarks in Section 3.4, describes a streaming algorithm that 
maintains an (9(l/e 1 / 2 )-size coreset at an amortized cost of (^((l/e 1 / 2 ) log 2 (l/e)) 
per insertion. For a fixed q, we go through the coreset (after computing its 
convex hull, if necessary), and determine the narrowest strip covering it and 
satisfying our angle constraints. The maximum of that and dist(g, £) gives the 
width of the minimum-width V-shape whose boundary passes through q. 3 The 
amortized cost per point is dominated by the (^((l/e 1 / 2 ) log 2 (l/e)) cost of 
insertion. Together with presorting points by distance from £, the total cost 
is then 0(nlogn + (n/e 1 / 2 ) log 2 (l/e)). □ 

Combining Lemmas 5 and 6 yields the procedure claimed in Lemma 4 and 
thereby completes our description of the (1 4- ^-approximation algorithm. 

7 Concluding remarks 

As mentioned in the introduction, this work was inspired by research on curve 
fitting, in the situations where a curve takes a sharp turn. Besides the exact 

3 More precisely, q lies on the boundary of Si and may not even appear on the boundary 
of V. However, as before, all V-shapes we examine are valid and cover P, and the 
desired approximating V-shape is among them, which is sufficient. 
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and approximate versions of the problem studied above, it would be natural 
to investigate a variant that can handle a small number of outliers. A natural 
"peeling" approach to the problem would be to eliminate the points defining the 
optimal V-shape found by our exact algorithm and trying again. However, it is 
easy to construct an example of a point set in which removing a single point not 
appearing on the boundary of the minimum- width covering V-shape significantly 
reduces the width of the optimum V-shape. 

Are there natural assumptions (perhaps in the spirit of "realistic input 
models" [9] or in the form of requiring reasonable sampling density) that would 
be relevant for the curve-fitting problem, and that would make finding the 
minimum-width covering V-shape easier? 

Returning to the problem studied in the paper, is it possible to find an 
exact minimum-width covering V-shape in subquadratic time? Is the problem 
3sUM-hard? 

Is it possible to speed up the approximation algorithm, improving the depen- 
dence of its running time on e? Is time 0(n + /(j)) achievable? 

Finally, we would like to point out that there are other "reasonable" definitions 
for a V-shape, if the goal is to approximate a sharp turn of a curve: One can 
imagine defining a V-shape as the Minkowski sum of a disk with the union of 
two rays emanating from a common point as in [17]. The width of the V-shape 
would be the diameter of the disk. Can the exact algorithm from [17] be sped 
up? Is there a faster approximation algorithm? Is this version of the problem 
better suited for curve fitting? 
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